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I.  INTRODUCTION 


The  ability  to  predict  structural  response  in  the  vicinity  of  a  large  and 
sudden  energy  source  is  dependent  upon  definition  of  the  loads,  the  behavioral 
characteristics  of  the  geological  material  surrounding  the  structure,  the 
mechanism  for  transferring  loads  from  the  geological  material  to  the  struc¬ 
ture,  and  finally  the  characteristics  of  the  structure.  The  basic  objective 
of  this  research  project  is  to  provide  an  improved  understanding  of  load 
transferral  from  one  medium  to  another  or,  more  specifically,  the  modeling  of 
structure-media  (concrete-soil)  interactions.  However,  to  address  this  prob¬ 
lem  it  is  essential  that  good  constitutive  models  be  available  for  both  soil 
and  concrete.  The  formulation  of  such  a  model  is  discussed  in  Section  II,  and 
the  research  efforts  that  are  needed  to  improve  the  model  are  described. 

The  plasticity  formulation  used  for  this  study  is  amenable  to  the 
incorporation  of  rate  effects  in  several  ways.  One  method  is  described  in 
Section  III,  and  theoretical  results  for  a  loose  sand  are  presented.  The 
model  shows  a  smooth  behavior  with  respect  to  the  two  parameters  that  are  used 
to  incorporate  rate  effects.  The  approach  is  therefore  believed  to  be 
adequate  for  now,  especially  in  light  of  the  paucity  of  experimental  data. 

This  plasticity  formulation  involves  a  parameter  that  can  be  interpreted 
as  the  shear  capacity  of  a  material  under  zero  mean  pressure.  With  a  slight 
modification,  this  parameter  can  be  used  to  provide  tensile  cutoff.  As  indi¬ 
cated  in  Section  IV,  the  algorithm  is  numerically  stable  for  such  an  applica¬ 
tion.  When  this  parameter  is  allowed  to  drop  to  zero,  the  incorporation  of 
fracture  into  the  model  becomes  a  realistic  possiblity. 

Concrete  and  many  geological  materials  exhibit  softening,  a  phenomenon 
that  causes  severe  numerical  problems  because  it  represents  a  state  in  which 
the  stiffness  of  a  system  ceases  to  be  positive  definite.  If  models  such  as 
the  one  described  in  this  report  are  to  be  used  in  large-scale  computer  codes, 
robust  and  efficient  algorithms  must  be  available  to  handle  both  dynamic  and 
static  loading  conditions.  Static  loading  conditions  are  especially  diffi¬ 
cult;  this  problem  is  discussed  briefly  in  Section  V. 
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ine  report  concludes  with  a  Drier  summary  of  the  resea 
associated  with  this  project.  Also  discussed  in  Section  VI 
areas  on  which  this  research  will  focus  as  the  work  continu 


II.  PLASTICITY  MODEL 


The  current  version  of  the  plasticity  model  is  described  in  Appendix  A. 

Of  particular  note  is  the  use  of  a  new  third  invariant  that  allows  the  limit 
surface  and  the  flow  surface  to  be  described  adequately  with  straight  lines. 
This  formulation  automatically  leads  to  a  good  representation  of  limit 
stresses  with  smooth  hardening  and  softening  features.  The  model  has  also 
been  used  to  provide  comparisons  for  a  variety  of  paths  and  materials  other 
than  those  shown  in  Appendix  A.  Detailed  comparisons  with  a  weak  concrete 
will  be  presented  in  a  future  report. 

One  of  the  fundamental  assumptions  implicit  in  the  theory  is  that  at 
the  limit  state,  some  measure  of  inelastic  strain  will  assume  the  same  value 
for  all  stress  paths.  The  current  measure  involves  a  path  length  in  the  space 
of  the  first  and  third  invariants  of  inelastic  strain.  It  is  found  that,  for 
example,  predictions  of  strain  at  the  limit  state  for  uniaxial  stress  and  pure 
shear  are  off  by  a  factor  of  2,  whereas  for  other  paths  the  predicted  strains 
are  within  5  percent  of  the  experimental  values.  These  results  indicate  the 
need  for  a  better  measure  of  inelastic  strain,  or  damage,  and  some  effort  will 
be  expended  on  this  task.  There  is  a  distinct  possibility  that  the  use  of  a 
single  measure  is  too  simplistic  and  that  two  or  more  parameters  will  have  to 
be  used. 

The  use  of  the  third  invariant  without  the  second  invariant  fixes  the 
shape  of  the  flow  surface  in  the  pi-plane.  The  flow  law  is  associated  in  this 
space,  and  for  the  majority  of  paths  considered,  the  assumption  is  good. 
However,  for  certain  proportional  loading  paths  it  appears  that  a  nonassoci- 
ated  law  may  be  required.  An  attempt  will  be  made  to  identify  closely  the 
conditions  under  which  the  associativity  assumption  must  be  relaxed. 

Little  thought  has  been  given  to  the  formulation  of  softening  within 
the  model.  However,  because  softening  can  be  controlled  through  the  relation 
involving  the  slope  of  the  flow  surface  and  the  inelastic  strain  measure, 
alternate  relations  for  simulating  observed  softening  phenomena  will  be 
considered. 


III.  RATE  EFFECTS 


The  formulation  given  in  Appendix  A  is  rate- independent.  A  conventional 
method  of  incorporating  rate  effects  is  to  use  an  overstress  approach  in  which 
the  rate  is  related  to  the  amount  by  which  a  measure  of  stress  exceeds  the 
same  measure  under  rate-independent  conditions.  Thus,  the  strain  rate  is 
governed  by  the  measure  of  stress.  An  alternate  approach,  which  is  the  one 
adopted  here,  specifies  the  yield  surface  in  terms  of  a  measure  of  inelastic 
strain  and  its  rate.  Here,  the  stress  is  governed  by  the  measure  of  strain 
that  is  used.  The  choice  of  one  approach  over  the  other  is  largely  a  matter 
of  preference  because  the  experimental  data  that  would  provide  the  basis  for  a 
selection  do  not  exist. 

The  slope  parameter,  y,  is  a  function  of  inelastic  measure  of  strain,  51. 
A  viscoplastic  theory  is  obtained  if  y  Is  replaced  by  a  parameter,  y  ,  that 
depends  on  strain  rate  as  well.  It  is  assumed  that  the  following  separable 
form  is  adequate: 

Yr  *  y  [l  -  A  tanh  (e/eo)  ]  (1) 

•  • 

In  which  A  is  an  amplification  factor,  e  is  an  inelastic  strain  rate,  and 
denotes  a  reference  strain  rate.  Sample  runs  were  performed  on  a  loose  sand 
for  which  experimental  data  on  the  static  case  are  available  (Ref.  1).  Fits 
to  these  data  are  shown  in  Figure  1  for  the  cases  of  tri axial  compression  and 
tri axial  extension.  These  tests  were  performed  in  a  cubical  tri axial  appara¬ 
tus  at  constant  confining  stresses  (a3  *  5882  Pa).  In  triaxial  compression 
the  intermediate  principal  stress,  ®2»  is  held  equal  to  the  confining  stress 
(o2  *  o3),  whereas  in  the  triaxial  extension  test  the  intermediate  and  major 
principal  stresses  in  compression  are  equal  («2  *  °i)*  These  tests  represent 
the  extreme  values  of  the  intermediate  principal  stress.  The  static  compari¬ 
son  is  good  except  in  the  case  of  the  volumetric  strain  associated  with  triax¬ 
ial  extension. 

•  • 

A  reference  strain  rate  of  1.0  was  chosen  for  Sq,  and  i  was  arbitrarily 
fixed  at  a  value  of  0.75.  The  effect  of  changing  A  was  evaluated  numerically. 
The  resu  *  »re  c  ,«n  In  Figure  1,  which  illustrates  that  the  response  changes 
smoothly  w'  this  amplification  factor  so  that  fitting  experimental  data 
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(a)  Triaxlal  compression. 

Figure  1.  Theoretical  predictions  for  a  range  of  rate  effects  and  static 
experimer  al  results  (Ref.  1)  for  loose  sand. 

8 


al  stress,  kg/cm 


should  be  relatively  straightforward.  The  present  formulation  indicates 
significantly  more  dilatation  with  a  strain  rate  effect.  Whether  or  not  this 
indication  is  realistic  remains  to  be  seen. 

The  only  data  that  involve  strain  rates  for  soils  in  the  range  of  inter¬ 
est  to  this  study  are  the  uniaxial  strain  data  for  a  given  stress  rate  (Refs. 
2  and  3).  These  data  will  be  used  to  evaluate  parameters  A  and  eo  as  an 
application  for  the  model  during  the  next  phase  of  this  research.  More 
effort  will  be  expended  on  the  formulation  of  a  suitable  experimental  program 
that  will  provide  the  data  needed  for  both  concrete  and  geological  materials. 


IV.  LIMITATIONS  ON  TENSILE  STRESS 


The  existing  formulation  is  limited  to  the  case  of  positive  mean  pressure 
because  not  all  of  the  implications  associated  with  the  new  third  invariant 
for  tensile  stresses  have  been  investigated.  However,  as  a  preliminary  inves¬ 
tigation  to  determine  whether  the  model  could  handle  tensile  cutoff,  a  flow 
rule  with  y  *  i*  4  0  or  <p  =  -L  -  o$  for  P  <  0  was  assumed.  The  strains  for 
uniaxial  tension  are  shown  in  Figure  2,  where  the  positive  strain  and  stress 
axes  denote  compressive  behavior.  As  the  uniaxial  stress  is  increased,  the 
corresponding  strain  also  increases  and  the  lateral  strains  indicate  contrac¬ 
tion  in  the  classical  Poisson  effect.  When  the  cutoff  stress,  o$J  is 
reached,  the  plasticity  model  automatically  yields  increasing  strain  with  no 
further  increase  in  stress,  while  the  lateral  strains  remain  fixed.  The 
uniaxial  strair  behavior  is  in  accordance  with  cracking  phenomena,  but  the 
stress  and  lateral  components  of  strain  should  drop  abruptly  to  zero. 

Predicted  stress-strain  behavior  for  a  biaxial  tensile  path  with 
a2  =  0.5  oi  is  shown  in  Figure  3.  Elastic  behavior  is  observed  until  the  state 
L  *  -<*s  (which  corresponds  to  Oj,  *  and  o2  =  0.5os)  is  reached.  For 
additional  straining  along  the  xraxis,  the  stresses  remain  constant  while  the 
lateral  strains  exhibit  a  slight  amount  of  extension.  This  behavior  is  not  a 
good  physical  representation,  but  it  may  be  a  reasonable  approximation  for 
postcracking  behavior. 

A  shear  tensile  path  defined  as  °i  *  -o2,  with  a3  *  0,  is  plotted  in 
Figure  4.  The  material  behaves  elastically  until  L  =  -o$  or  ^  =  as  and 
°2  =  after  which  the  stresses  remain  constant,  the  tensile  strain  con¬ 
tinues  to  increase,  and  the  other  strains  do  not  change.  This  is  a  realistic 
representation  of  physical  behavior. 

These  examples  Indicate  that  with  a  minor  modification  to  the  flow  sur¬ 
face,  tensile  cutoff  can  be  represented  within  the  existing  procedure.  The 
algorithm  remains  stable  as  one  strain  continues  to  increase  in  a  monotonic 
fashion.  However,  the  approach  is  too  simplistic  to  represent  cracking 
because  most  data  seem  to  Indicate  that  cracking  is  Initiated  when  the  maximum 
principal  stress,  rather  than  the  mean  stress,  reaches  a  critical  value. 


VSP. 


The  prediction  of  reasonable  strains  in  the  cutoff  region  can  be  attrib¬ 
uted  to  the  fact  that  the  stress,  c$,  remains  finite.  On  the  other  hand,  a 
more  realistic  representation  for  stresses  is  obtained  if  o$  drops  to  zero; 
thus,  the  requirements  for  stress  and  strain  appear  to  set  up  conflicting 
demands  on  the  model.  It  is  believed  that  the  issue  can  be  handled  neatly  by 
noting  that  with  a  nonassociated  flow  rule,  a  considerable  degree  of  flexibil¬ 
ity  is  available.  Thus  it  is  proposed  that  the  flow  condition  and  the  poten¬ 
tial  function  be  modified  to  the  following  forms: 


-L(°s)  +  YP  -  °s 

*  *  *  * 
*  3  -L  °s  +  1  P  - 


(2) 


in  which  o*  is  set  equal  to  for  the  beginning  state.  Then  if  a  crack  is 

predicted  according  to  some  criterion,  a  will  be  reduced  to  zero,  which  will 

^  * 

yield  an  appropriate  stress  field,  while  o$  will  be  increased,  which  will 
yield  an  appropriate  displacement  field.  This  condition  can  be  accomplished 
with  very  little  modification  to  the  existing  algorithm. 


Even  this  approach  will  not  yield  the  directional  properties  associated 

with  cracking.  To  obtain  these  properties,  it  is  proposed  that  the  more  gen 

* 

eral  form  of  <j»  and  $  ,  in  which  the  tensor  o  is  used  instead  of  oK  be 

'"5 

investigated.  The  questions  to  be  answered  are  (1)  whether  the  approach  is 
feasible,  and  (2)  if  It  is  feasible,  whether  the  procedure  offers  any  advan¬ 
tages  over  the  use  of  a  cracking  algorithm. 
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V.  NUMERICAL  SOLUTION  PROCEDURES 


A  constitutive  model  that  yields  correct  response  features,  including 
softening  for  large  strains,  provides  a  means  of  using  less  conservative 
designs  based  on  the  actual  energy-absorbing  capacity  of  frictional  materials. 
These  response  features  involve  stress  paths  that  go  beyond  the  limit  state. 

To  predict  behavior  into  this  regime,  a  robust  solution  algorithm  is  required. 
Several  exist  for  dynamic  problems.  The  reason  is  that  the  mass  matrix 
remains  positive  definite;  consequently,  it  is  not  difficult  to  obtain  numeri¬ 
cal  approximations  to  the  acceleration  from  the  equation  of  motion.  However, 
the  situation  is  quite  different  for  softening  behavior  with  problems  that 
involve  static  loads  or  those  that  reflect  transition  from  a  dynamic  to  a 
static  state.  An  example  of  the  latter  case  is  a  structural  member  subjected 
to  an  impulse  sufficiently  large  to  force  the  response  into  the  softening 
regime.  For  either  case  the  tangent  stiffness  matrix  ceases  to  be  positive 
definite,  and  conventional  solution  procedures  are  no  longer  applicable. 
Alternate  algorithms  are  available,  but  convergence  and  efficiency  are  matters 
of  concern. 

One  approach  that  has  considerable  appeal  is  dynamic  relaxation  (Ref. 

4),  which  is  an  explicit  iterative  method  for  obtaining  static  solutions.  It 
is  based  on  the  fact  that  the  static  solution  is  the  steady-state  part  of  the 
transient  response  for  a  temporal-step  load,  and  it  uses  the  familiar  dynamics 
equations  with  damping.  The  amount  of  damping  is  based  on  the  critical  damp¬ 
ing  factor  associated  with  the  fundamental  frequency,  which  in  turn  is  related 
to  the  lowest  eigenvalue  of  the  tangent  stiffness  matrix.  If  this  lowest 
eigenvalue  is  zero  or  negative  (strain  softening),  the  damping  factor  is  zero. 
For  this  reason  the  numerical  solution  converges  very  slowly  under  these  con¬ 
ditions.  For  more  conventional  strain  hardening,  the  procedure  yields  a  solu¬ 
tion  with  an  efficiency  similar  to  that  obtained  by  using  conventional  matrix 
equation  solvers. 

Padovan  (Ref.  5)  has  recently  developed  a  procedure  that  Involves  the 
use  of  limitations  on  a  norm  of  the  displacement  and  applied  load  vectors. 
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With  this  measure  he  is  able  to  track  the  local  nature  of  the  response  (hard¬ 
ening  or  softening)  and  place  suitable  limitations  on  increments  in  either 
force  or  displacement  to  ensure  accuracy  and  convergence.  The  procedure  has 
been  applied  to  a  variety  of  problems  (Ref.  6),  but  it  is  not  clear  whether 
the  algorithm  is  suitable  for  use  in  a  large-scale  computer  code.  Neverthe¬ 
less,  the  approach  appears  to  be  a  distinct  improvement  over  existing  methods, 
and  it  must  be  considered  seriously. 

Quasi-Newton  methods,  of  which  BFGS  is  a  special  case,  are  also  being 
applied  to  problems  of  this  type  (Refs.  7  and  8).  These  methods  involve  var¬ 
ious  techniques  for  numerically  updating  the  Jacobian,  not  necessarily  at 
every  step,  to  provide  an  appropriate  Newton  step  for  minimizing  a  functional. 
Both  local  and  global  strategies  are  used  to  ensure  convergence,  and  various 
steps  are  taken  to  preserve  symmetry  and  sparsity.  Also,  the  step  size  and 
options  for  backtracking  are  chosen  for  optimal  efficiency.  Although  these 
concepts  are  still  in  a  developmental  mode  as  far  as  engineering  computations 
are  concerned,  it  is  believed  that  many  of  them  will  prove  to  be  directly 
applicable  to  engineering  problems  that  exhibit  softening. 

An  avenue  that  seems  to  hold  considerable  promise,  but  which  apparently 
has  not  been  explored,  is  the  combination  of  some  of  the  concepts  of  Padovan 
and  the  quasi-Newton  procedures  with  dynamic  relaxation.  Because  of  its  com¬ 
putational  simplicity  for  large-scale  problems  and  because  basic  equations  of 
dynamics  are  used,  the  dynamic  relaxation  procedure  has  considerable  appeal  to 
engineers.  Furthermore,  it  forms  a  natural  procedure  for  making  the  transi¬ 
tion  from  dynamic  to  static  problems,  even  in  a  softening  regime.  Because 
slow  convergence  is  a  critical  matter  for  large  problems,  this  would  seem  to 
be  a  particularly  productive  research  area  in  light  of  the  need  for  reliable 
and  efficient  numerical  algorithms. 
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VI.  RESEARCH  TOPICS 


The  primary  aim  of  this  project  to  date  has  been  to  determine  the  feasi¬ 
bility  of  using  a  simple  plasticity  model  for  frictional  materials.  It  is 
believed  that  enough  evidence  has  been  accumulated  to  demonstrate  conclusively 
that  features  such  as  strain  hardening,  shear  enhanced  compaction,  dilatation, 
strain  softening,  and  limit  states  evolve  naturally  from  the  model.  An  asso¬ 
ciated  numerical  algorithm  is  also  functioning  and  appears  to  be  reliable  and 
fairly  efficient. 

The  primary  problem  associated  with  the  model  is  the  accuracy  of  the 
strain  predictions.  For  many  paths,  strains  are  within  ±10  percent  of  the 
experimental  values  for  all  stresses,  which  is  considered  more  than  adequate 
for  engineering  analyses  involving  these  materials.  On  the  other  hand,  pre¬ 
dicted  strains  at  the  limit  state  for  some  paths  differ  from  experimental 
values  by  a  factor  of  2.  This  discrepancy  is  not  acceptable  and  can  be 
attributed  to  the  choice  of  a  particular  inelastic  strain  invariant.  One 
objective  of  the  continuing  effort  will  be  to  examine  closely  the  available 
experimental  data  in  order  to  determine  whether  a  better  invariant  can  be 
defined.  The  objective  is  to  define  an  invariant  in  such  a  way  that  the  value 
of  the  invariant  at  the  limit  state  is  the  same  for  all  paths. 

Preliminary  calculations  involving  strain  rate  effects  have  been  per¬ 
formed.  This  exploration  will  continue,  as  will  the  accumulation  of  experi¬ 
mental  data— especially  those  data  for  multiaxial  states  of  stress  that  show 
general  trends.  The  objective  here  is  to  attempt  to  identify  the  functional 
form  involving  strain  rates  for  flow  surfaces  that  is  best  suited  for  matching 
experimental  data.  To  enhance  the  experimental  data  base,  an  experimental 
program  for  concrete  will  be  defined,  and  then  a  testing  program  will  be 
Implemented. 

A  potential  cracking  model  that  fits  in  with  the  plasticity  formulation 
has  been  identified.  If  this  approach  should  prove  feasible,  one  constitutive 
relation  and  algorithm  would  handle  both  the  compressive  and  the  tensile 
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behaviors  of  frictional  materials.  The  computational  advantages  would  be 
considerable.  Furthermore,  the  formulation  would  provide  a  natural  framework 
for  handling  soil-structure  interfaces.  The  anisotropic  feature  that  is  auto¬ 
matically  present  with  the  existence  of  a  crack  could  be  exploited  to  test 
postulates  on  how  shear  and  normal  stresses  interact  at  an  interface.  Again, 
carefully  controlled  laboratory  data  will  be  sought  to  provide  guidelines  on 
suitable  approaches. 

The  representation  of  strain  softening  has  generally  been  overlooked  in 
determining  the  energy- absorption  capabilities  of  frictional  materials. 

Because  most  numerical  algorithms  cannot  cope  with  softening,  an  artificial 
limit  stress  is  often  used  instead  of  the  actual  one.  With  such  an  approach, 
there  is  no  way  to  accurately  ascertain  whether  an  approach  is  conservative, 
and  a  great  deal  of  engineering  judgment  is  required.  Strain  softening  is 
controlled  through  the  use  of  an  inelastic  strain  invariant  in  the  plasticity 
model;  hence,  a  method  of  matching  data  is  available.  This  aspect  of  the 
model  will  be  explored  to  a  limited  extent  to  determine  exactly  what  data  are 
required  and  why. 

* 

Ultimately,  any  good  constitutive  relation  will  be  used  in  large-scale 
computer  codes.  Because  the  plasticity  model  can  accomodate  strain  softening, 
which  is  a  localized  feature,  it  becomes  feasible  on  a  modeling  basis  to  cap¬ 
ture  softening  on  a  global  scale.  One  example  is  the  response  to  failure  of  a 
beam  subjected  to  a  lateral  load  that  is  displacement-controlled.  No  cur¬ 
rently  available  equation  solver  handles  such  problems  routinely.  The  tangent 
stiffness  method  becomes  singular,  and  the  dynamic  relaxation  procedure  is  too 
inefficient  for  application  to  engineering  problems.  A  numerical  algorithm 
that  could  handle  this  problem,  as  well  as  dynamic  problems,  is  needed.  Moni¬ 
toring  of  the  literature  will  be  continued  so  that  new  developments  can  be 
investigated  as  the  need  arises. 

Although  several  specific  topics  have  been  addressed  in  this  section, 
they  are  all  unified  around  the  single  theme  of  developing  accurate  models  for 
simulating  soil -structure  Interaction.  This  focus  will  be  maintained,  with 
special  emphasis  on  strain  invariants,  rate  effects,  and  cracking. 
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ABSTRACT 


A  nonassoci ated  plasticity  theory  is  given  for  frictional  materials  in 
terms  of  the  first  and  third  invariants  of  stress  and  inelastic  strain.  Fea¬ 
tures  such  as  strain  hardening,  strain  softening,  dilatation,  and  compaction 
are  exhibited.  Limit  states  and  the  Ko-condition  are  interpreted  in  a  natural 
manner  with  the  theory.  Comparisons  of  predictions  with  experimental  data  for 
limit  states  and  deformation  paths  are  given  for  several  materials. 


♦This  appendix  is  a  reproduction  of  a  professional  paper  accepted  for  pub¬ 
lication  in  Journal  of  Structural  Mechanics.  Thus,  it  is  a  self-contained 
document  with  its  own  internally  consistent  numbering  system  for  equations, 
references,  and  figures  and  with  a  format  prescribed  by  the  publisher. 
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I.  INTRODUCTION 


The  conventional  theory  of  plasticity  has  a  long  history  of  successful 
use  In  applications  Involving  metals.  The  theory  Is  well-founded  on  a  phys¬ 
ical  and  a  mathematical  basis  with  numerous  applications  to  engineering 
problems  through  the  use  of  computer  programs.  This  wealth  of  experience 
and  familiarity  has  led  to  numerous  attempts  to  apply  the  theory  to  fric¬ 
tional  materials  such  as  concrete,  rock,  soils.  Ice,  and  snow.  Although 
there  Is  no  fundamental  reason  why  the  theory  of  plasticity  should  be  appro¬ 
priate  for  these  materials,  there  are  now  several  formulations  that  provide 
correct  behavioral  characteristics.  Unfortunately  many  of  these  relations 
are  so  complicated  that  their  use  Is  limited  to  the  simplest  engineering 
applications.  This  paper  Is  the  result  of  an  attempt  to  formulate  a  model 
that  captures  essential  response  characteristics  without  recourse  to  a  large 
number  of  parameters  or  special  conditions.  The  Introduction  of  a  new  third 
Invariant  leads  to  a  plasticity  theory  that  appears  to  be  natural  for  fric¬ 
tional  materials  In  the  same  sense  that  the  von  Mlses  formulation  works  so 
well  for  a  wide  class  of  metals. 

Drucker  and  Prager  [1]  generalized  the  Mohr-Coulomb  hypothesis  for  slip 
to  a  conical  yield  function  consisting  of  a  linear  combination  of  the  first 
and  second  Invariants.  This  function  was  then  used  to  provide  upper  and 
lower  bounds  for  the  critical  height  of  a  vertical  bank  with  certain  Impli¬ 
cations  concerning  sliding  surfaces  and  volume  expansion  (dilatation).  This 
theory  of  perfect  plasticity  Is  too  Idealized  to  adequately  represent  the 
behavior  of  soils  so  the  concept  of  work-hardening  was  applied.  To  account 
for  volumetric  compaction  under  hydrostatic  loads  and  to  control  the  amount 


of  dilatation,  a  cap  was  placed  on  the  open  end  of  the  cone  by  Drucker,  Gib¬ 
son,  and  Henkel  [2].  Modifications  to  this  model  have  been  made  over  the 
years  and  the  fact  that  the  approach  Is  still  In  use  [3]  Is  testimony  to  the 
Insight  provided  by  these  papers. 

Inevitably,  as  more  experimental  data  became  available,  the  limitations 
of  the  assumed  shape  for  the  yield  surface  became  apparent.  In  the  pi-plane 
the  circular  shape  of  von  Mlses  must  be  replaced  by  a  more  triangular  shape 
as  expressed  by  the  parametric  form  of  Wiliam  and  Warnke  [4],  This  general 
shape  Is  also  obtained  If  the  third  Invariant  of  stress  Is  used,  either 
directly  as  suggested  by  Lade  [5]  and  Vermeer  [6],  or  Indirectly  through  the 
use  of  the  Lode  angle  as  exemplified  by  the  work  of  Mroz,  Norris,  and  Zlen- 
klewlcz  [7]  and  Hsleh,  Ting,  and  Chen  [8].  More  sophisticated  versions 
Involving  multiple  surfaces  have  been  proposed  [9-12]  primarily  to  handle 
cyclic  loading  In  a  smooth  manner  consistent  with  experimental  data.  Of 
course,  numerous  rate  type  models  have  also  been  suggested  for  frictional 
materials  [e.g.,  13-15],  but  the  requirement  of  a  large  number  of  material 
parameters  has  prevented  their  adoption  to  any  significant  scale. 

The  majority  of  these  latter  models  are  similar  In  that  they  Involve 
all  three  stress  Invariants  and  that  results  are  presented  In  a  plane 
Involving  the  first  and  second  stress  Invariants.  To  a  certain  extent  this 
has  produced  a  degree  of  complexity  that  Is  not  present  In  the  original  for¬ 
mulations  of  References  1  and  2.  Although  computers  provide  the  luxury  of 
accommodating  several  parameters,  limits  on  storage  are  quickly  reached  for 
many  engineering  applications  so  that  there  Is  still  a  great  need  for  formu¬ 
lations  that  capture  essential  response  features  but  remain  simple  In  form. 
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Lade  and  Duncan  [16]  observed  that  the  limit  state  for  a  coheslonless 
soil  can  be  represented  without  the  use  of  the  second  Invariant  and  Lade 
[17]  showed  recently  that  concrete  also  falls  In  this  category.  This  sug¬ 
gests  that  a  natural  plasticity  formulation  for  frictional  materials  should 
Involve  only  the  first  and  third  Invariants  rather  than  all  three.  Lade  [5] 
has  developed  this  Idea  to  a  certain  extent  but  he  transforms  all  of  his 
results  Into  conventional  spaces,  which  Is  done,  naturally  enough,  to  facil¬ 
itate  communication  and  comparison  with  previous  work.  However,  It  Is 
believed  that  In  attempting  to  use  conventional  concepts,  the  natural  ele¬ 
gance  of  the  theory  Is  obscured.  Here  an  attempt  Is  made  to  present  a  the¬ 
ory  In  as  simple  a  mathematical  form  as  possible  with  some  comparison  with 
data  to  Illustrate  that  qualitatively  correct  response  features  are  obtained. 


2.  LIMIT  SURFACES  AND  STRESS  PATHS 

Limit  surfaces  for  concrete,  rock,  soils.  Ice,  and  snow  depend  on  tfye 
first  Invariant  of  stress,  or  mean  pressure,  and  are  known  as  'rlctlonaf 
materials  In  contradistinction  to  most  metals.  These  surfaces  also  depend 
on  a  measure  of  shear  which  can  be  expressed  through  the  second  and  third 
Invariants  of  stress.  The  conventional  approach  has  been  to  express  a  limit 
surface  In  terms  of  the  mean  pressure  and  the  second  Invariant  of  the  stress 
deviator 
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with  modifications  provided  by  Incorporating  a  third  Invariant  as  dictated 
by  experimental  observations.  Lade  and  Duncan  [16]  and  Lade  [17]  have  shown 
that  If  a  particular  form  of  the  third  Invariant  Is  used,  the  second  Invari¬ 
ant  Is  not  required  In  the  expression  for  a  limit  surface  for  at  least  par¬ 
ticular  classes  of  soils  and  concrete.  A  modification  of  their  approach  Is 
proposed  to  utilize  this  feature,  which  may  hold  for  all  frictional 
materials . 

With  the  use  of  a  shift  In  stress,  o,  a  convenient  definition  for  a 
third  Invariant  Is 


L  ■  -  [det  (a  -  a  )  +  det  (or  )]  x/ *  (2) 

~  ~s  ~s 

In  which  the  second  term  In  the  brackets  Is  Included  to  make  L  equal  to  zero 
for  a  state  of  zero  stress.  Analogous  to  the  square  root  that  Is  often  used 
with  the  second  Invariant,  a  cube  root  Is  Introduced  to  provide  the  dimen¬ 
sion  of  stress.  For  principal  stresses  ax,  oa,  a >  and  an  Isotropic 

assumption  for  the  shift  stress,  a  -  a  I,  the  first  and  third  Invariants 

~s  s  - 

reduce  to 


P  «  -  ”  (ffi  ♦  #»  f  ffi) 

L  -  -  -  a  )  (9a  -  o$  )  ( 9 1  -  os)  +  o*]x/* 


(3) 


A  limit  point  Is  the  point  on  a  prescribed  stress  or  strain  path  at 
which  the  state  of  stress  Is  stationary  with  respect  to  Increments  In 
strain.  If  such  points  for  a  variety  of  paths  are  plotted  In  stress  space, 
the  result  of  Interpolating  between  the  points  Is  a  limit  surface  (often 
called  a  failure  surface).  Limit  points  for  concrete,  granite,  marble, 
dense  and  loose  sands,  and  clay  are  shown  In  Figure  1.  A  reasonable  fit  to 
these  points  Is  postulated  to  be  the  line 
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(a)  Low-  and  high-strength  concrete. 


Figure  1.  Limit  curves  for  frictional  materials. 
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A  Grundlte  clay.  Ref.  21 
O  Loose  sand.  Ref.  16 
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Figure  1.  Concluded 


(4) 


described  with  two  parameters  consisting  of  a  limit  slope,  and  the 
Intercept  on  the  L-axIs  represented  by  the  shift  stress,  a^.  The  fits 
are  remarkably  good  considering  the  elementary  formulation.  The  scatter 
that  Is  present  Is  not  unusual  for  these  materials  and  Indicates  the  need 
for  probabilistic  theory  If  a  proper  engineering  analysis  Is  to  be  performed. 

The  use  of  a  limit  surface  In  the  P-L  space  suggests  that  this  might 
also  be  a  convenient  space  to  formulate  flow  surfaces  for  use  with  the  the¬ 
ory  of  plasticity.  However,  this  Is  an  unconventional  space  and  the  charac¬ 
ter  of  simple  stress  paths  Is  not  Intuitively  obvious.  Also,  there  Is  con¬ 
siderable  benefit  to  be  able  to  relate  back  to  other  formulations  and  to 
show  the  potential  advantages  of  working  In  this  space. 

In  principal  stress  space,  consider  a  transformation  to  stresses 
defined  by 

q*.  »  {ax  -  aa)/^2 

qa  ■  (2oi  -  «i  -  oa)/yf 6*  (5) 

qa  -  -  {ax  ♦  9a  ♦  o*)A/T 

The  qs-axls  defines  the  direction  of  hydrostatic  pressure  and  the  qa.  - 
qa  plane  Is  called  the  pi-plane.  Stress  paths  and  limit  surfaces  are  fre¬ 
quently  shown  In  the  pi-plane  and  such  a  plot  Is  Included  for  the  sake  of 
completeness . 

To  Illustrate  stress  paths  In  the  P-L  space,  let  a  denote  a  positive, 
monotonlcally  Increasing  parameter.  Then  a  series  of  conventional  load 
paths  are  defined  as  follows  for  P  >  0: 
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(a)  Hydrostatic  compression 
0i  m  a  a  m  0i  ■  -a 

P  -  a  L  -  [{a  a  )*-  crs*]  l/# 

qx  -  0  qa  -  0 

(b)  Uniaxial  compression 

<7a  *  -a  oi  ■  oi  ■  0 

_  .  a  .  x  /  a 

P  ■  a/3  L  -  (aos  ) 

qx  «  a/)/ 2  qa  »  a//6 

(c)  Trlaxlal  compression  (hydrostatic  compression  to  P  >  a« 
followed  by  uniaxial  compression) 

01  «  -ao  -  a  0i  ■  os  ■  -ao  ao  >  a 

2  i  x  / 

P  -  ao  ♦  a/3  L  ■  [(ao  t  a  ♦  0s)|ao  *  ^  ~  ffs  1 

qx  -  a//T  qa  -  a//T 

(d)  Trlaxlal  extension  (hydrostatic  compression  to  P  -  ao  fol 
lowed  by  uniaxial  release) 

<7 x  m  0j  m  -ao  o a  ■  -ao  +  a  ao  >  a 

P  -  ao  -  a/3  L  ■  [(ao  -  a  *  °  s)i  ao  +  a*  )  -  <7*  ]l/* 

qx  -  0  qa  -  2a/v^6 

(e)  Biaxial  compression 

ox  ■  a  a  *  —a  01  ■  0 

P  -  2 a/3  L  -  [(a  ♦  a)a  o  *  -  o*s ]a/* 

qx  •  0  qa  *  2a//6" 
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(f)  Pure  shear  (hydrostatic  compression  to  P  •  ao  followed  by 
shear ) 

ax  -  -ao  a  ai  •  -ao  -  a  o»  -  0 

„  .  .  ,  a  .  x  /  a 

P  -  0  L»  -  -(a  <r s ) 

qi  -  2y*  qa  -  0 

These  paths  are  shown  In  Figure  2  for  -  1  and  In  the  pi-plane  of 
Figure  3. 

The  representation  of  the  limit  surface  In  the  pi-plane  Is  shown  In 
Figure  3  for  various  values  of  mean  pressure.  For  low  values  of  P  the  char¬ 
acteristic  triangle  with  rounded  corners  Is  obtained,  and  for  large  P  the 
surface  becomes  circular,  which  Is  a  feature  that  Is  generally  observed  for 
frictional  materials. 

Lade  associates  o  with  limiting  tensile  behavior.  The  formulation 
s 

here  suggests  an  alternate  Interpretation.  Consider  the  pure  shear  loading 

path  with  ao  ■  0  and  a  -  a  .  Then  L  •  -o  and,  since  this  Is  a  point  on 

s  s 

the  limit  surface  (actually  by  design),  a s  can  be  Interpreted  as  the 
limit  stress  of  a  material  under  pure  shear  with  zero  mean  pressure. 

A  glance  at  the  stress  paths  shown  In  Figure  2  shows  that  the  paths 
fall  between  the  hydrostatic  path  and  the  limit  surface.  In  other  words, 
conventional  stress  paths  fall  within  a  relatively  narrow  band  In  this 
space,  which  Is  simply  a  reflection  of  the  similarity  of  response  curves 
observed  for  this  class  of  materials.  Certain  stress  paths,  such  as  trlax- 
lal  compression,  become  almost  tangent  to  the  limit  surface,  so  It  Is  not 
surprising  that  experimental  data  for  limit  states  show  considerable  scat¬ 
ter.  For  large  values  of  P,  slightly  different  values  for  yield  rel¬ 
atively  large  differences  In  L  and  these  differences  can  be  magnified  even 
further  If  the  limit  points  are  plotted  In  a  different  space. 


Figure  2.  Stress  paths  and  flow  surfaces  In  the  P-L  plane 


A  uniaxial  strain  path  Is  one  of  the  conventional  experimental  paths 
used  for  soils.  As  the  strain  Is  Increased,  the  lateral  stresses  approach  a 
factor,  Ko,  times  the  applied  stress  and  this  Is  called  the  “Ko  condi¬ 
tion.*  In  the  P-L  space  this  path  depends  on  the  material,  but  It  Is  rea¬ 
sonable  to  suppose  that  the  path  asymptotically  approaches  the  limit  sur¬ 
face.  For  an  applied  stress  01  -  -a,  and  If  a  »  o^,  then  In  the  Ko 
condition  era  -  o»  -  -Koa  and 


P  -  J  (1  ♦  2Ko) 

,  1/  a/* 

L  *  cl  Ko 


(6) 


so  that  at  the  limit  surface 


The  Implication  Is  that  If  Ko  Is  Known  for  a  particular  material,  an 
estimate  for  Is  available  and  vice  versa.  For  example.  If  Ko  - 
0.40  then  -y^  -  0.90.  Aiso,  the  limiting  case  of  Ko  «  1,  ^  »  1 
might  be  used  to  represent  liquefaction. 


3.  FLOW  SURFACES  AND  FLOW  RULES 

For  the  development  of  a  plasticity  relation,  the  Initial  flow  surface 
can  be  developed  from  the  point  of  deviation  from  linearity  on  stress-strain 
curves.  Then  this  yield  surface  can  be  assumed  to  evolve  out  to  a  limit 
surface  to  represent  the  strain  hardening  phase  of  deformation.  Strain 
softening  can  be  obtained  by  letting  the  flew  surface  collapse  In  a  con¬ 


trolled  manner . 


Since  the  limit  surface  Is  conveniently  represented  by  a  line  In  a 
space  Involving  the  third  Invariant,  a  flow  surface  of  a  similar  nature  can 
be  postulated  as  follows: 

♦  -  yP  -  L  -  a  (8) 

In  which  y  Is  a  strain-hardening  function  with  the  restriction  y  >  yL. 

A  surface  *  -  0  defines  the  flow  state  and  Is  shown  as  a  dotted  line  In 
Figure  2.  For  stress  states  above  this  line,  4>  <  0,  which  represents  an 
elastic  state  and  *  >  0  Is  not  permitted. 

A  general  hardening  description  Is  obtained  by  allowing  the  flow  sur¬ 
face  to  rotate  toward  the  limit  line,  l.e.,  to  let  y  decrease  with  some 
measure  of  Inelastic  strain  to  the  value  y  .  Then,  with  further  defor¬ 
mation,  the  flow  surface  rotates  back  and  this  represents  strain  softening. 
The  specific  character  of  a  material  Is  expressed  In  the  functional  form  for 
y .  Some  materials  exhibit  very  little  softening  In  which  case  y  would 
Increase  only  slightly  after  reaching  the  limit  state.  Extreme  softening  Is 
obtained  by  allowing  y  to  Increase  significantly  after  the  limit  state. 

Most  plasticity  formulations  for  frictional  materials  are  presented  In 
terms  of  the  first  and  second  Invariants  of  stress.  For  convenient  compari¬ 
sons,  the  flow  surface  of  Equation  (8)  Is  shown  In  Figure  4  for  qa  -  0  and 
various  values  of  y.  For  y  >  1  rather  conventional  elliptical  shapes 
are  automatically  obtained  together  with  a  general  hardening  behavior 

that  would  be  extremely  difficult  to  develop  In  the  P  plane.  This  Is 

a 

another  example  that  Indicates  that  the  use  of  a  i.lrd  Invariant  Is  natural 
for  frictional  materials.  For  large  P  the  slope  of  the  hydrostatic  path  is 
unity,  so  If  y  >  1  the  flow  surface  does  not  Intersect  the  hydrostatic 
path.  This  Is  reflected  as  a  curve  In  the  p  plane  that  does  not  Inter- 

a 


sect  the  P-axIs  and  Is  a  consequence  of  the  simplest  possible  assumption 
that  the  flow  surface  Is  a  straight  line  In  the  P-L  plane. 

One  of  the  fundamental  difficulties  with  a  plasticity  formulation  Is 
that  appropriate  measures  of  deformation  must  be  selected  for  use  In  a 
strain  hardening  law.  If  these  measures  are  chosen  to  be  the  Inelastic 
strain  tensor,  e^ ,  then  the  theory  Is  usually  too  complicated  to  be  of 
much  practical  use,  so  a  positive  definite  Invariant,  e\  Is  often  adopted. 

If  a  reference  value,  e\  can  be  associated  with  a  limit  state  for  all  paths, 

O 

a  hardening  law  that  provides  the  appropriate  characteristic  features  Is 
given  by 


Y  -  r# 


( Y*  -  Yl  )  sin  j 


n 


(9) 


where  -yo  Is  the  value  of  y  that  defines  the  Initial  yield  surface.  The 
exponent,  n,  merely  controls  the  rate  at  which  y  changes  with  e1  . 

In  metal  plasticity,  e1  Is  often  chosen  to  be  the  path  length  associ¬ 
ated  with  the  second  Invariant  of  the  Inelastic  strain  deviator.  Here,  both 
volumetric  and  shear  effects  are  Important,  so  It  Is  postulated  that  an 
appropriate  Invariant  Is  a  path  length  Involving  all  three  invariants: 


-1 

e  - 

J* 

I  de*  I  +  b2del  +  ba  del 

*  O 

where  s  Is  a  monotonlcal ly  Increasing  load  parameter,  b2  and  bs  are  material 
constants,  and 


de1  ■  tr  (de1) 

V 

de]  =  [tr  (de1d  deid)] 

i  .•  1/3 

de]  »  | detCde1 ) | 


(11) 
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represent  Inelastic  volumetric  and  shear  Increments  based  on  the  Inelastic 


strain  tensor  and  Its  deviator,  e^d.  Such  a  formulation  has  proved  to  be 
useful,  at  least  for  particular  cases  of  sand  and  concrete,  but  It  must  be 
considered  as  a  preliminary  form  for  Illustrative  purposes. 

If  an  associated  flow  rule  Is  adopted,  then 

1  a* 

d  e  -  dX^J  (12, 

for  a  monotonlcally  Increasing  parameter  X  with  the  Increment  dX  obtained 
by  satisfying  the  consistency  condition,  $  C.  However,  In  agreement 
with  many  other  Investigators,  such  a  flow  rule  yields  an  excessive  amount 
of  dilatation,  so  the  use  of  a  potential  function,  $*,  Instead  of  $  Is 
required  In  Equation  (12).  After  an  Investigation  of  several  stress  paths, 
a  modification  to  y  was  deemed  sufficient  and  hence 

♦*  -  y*  P  -  l  -  a  (13) 

was  chosen  with  the  nonassoclated  parameter  y*  given  by  the  simple  additive 
relation 

Y*  =  Y  +  bi  (14) 

where  the  constant  bi  depends  on  the  material.  The  Implication 

of  Equation  (13)  Is  that  the  flow  rule  Is  associated  with  respect  to  L  but 

not  with  respect  to  P. 

In  terms  of  principal  stresses,  the  flow  law  yields  the  following 
expression  for  the  Inelastic  volumetric  strain  Increment: 

i  d\  r 

d  %  “  3Lal(<ri  "  °s)(aa  '  ffs  >  +  (<Ta  -  ffs  )(ffa  ‘  °s  > 

♦  (0*  -  <Ts  )  (01  -  0S  )  -  3La  Y*]  (15) 
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The  value  of  y*  that  gives  the  transition  point  between  dilatation  and 
compaction  Is  found  by  satisfying  the  condition  d  e^  -  0.  For  example,  under 
uniaxial  compression  with  the  stress  parameter  a  defined  previously,  this 
value  Is 


(16) 


(17) 


which  provides  some  general  guidance  as  to  how  the  function  y*  should 
behave  If  experimental  data  are  available  for  the  path.  For  hydrostatic 
compression,  y*  •  1 ,  and  since  compaction  Is  always  observed  for  such  a 
state,  the  restriction  y*  >  1  Is  Imposed  In  Equation  (14). 

Additional  Insight  Into  the  required  behavior  for  y*  Is  obtained  by 
considering  uniaxial  strain.  Since  compaction  Is  obtained  for  this  path, 
y*  must  be  greater  than  y*.  However,  at  the  Ko  condition,  the  slope 
of  the  mean  pressure-volumetric  strain  curve  Is  very  steep  which  Implies 
very  small  Increments  In  Inelastic  deformation.  Since  the  uniaxial  strain 
component  equals  the  volumetric  strain  for  this  case,  y*  must  be  greater 
than  but  close  to  y*.  The  use  of  the  Ko  stress  state  and  Equation  (15) 
yields 


y 


* 

t 


2  »  Ko 
3Kol/# 


(18) 


which  gives  y*  ■  1.09  for  Ko  -  0.4.  This  Is  the  type  of  Information  that 
Is  useful  for  determining  the  parameters  bi,  ba  and  ba. 
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4.  COMPARISON  WITH  EXPERIMENTAL  DATA 


To  demonstrate  the  flexibility  of  the  model  two  materials  with  consid¬ 
erably  different  response  characteristics  were  chosen.  The  first  Is  a  high- 
strength  concrete  that  exhibits  relatively  little  compaction  but  a  consider¬ 
able  amount  of  dilatation.  The  second  Is  a  material  similar  to  sand  but 
with  a  large  amount  of  porosity  so  that  severe  compaction  but  little  dilata¬ 
tion  Is  exhibited.  For  both  cases,  data  were  obtained  from  hydrostatic  and 
trlaxlal  compression  tests.  The  plots  of  the  trlaxlal  results  show  stress 
and  strain  measured  relative  to  a o  and  eo  which  denote  the  maximum  value 
for  stress  and  strain,  respectively,  obtained  In  the  Initial  hydrostatic 
part  of  the  path. 

Results  for  the  high-strength  concrete  tested  by  Green  and  Swanson  [18] 
are  shown  In  Figure  5  with  the  following  material  parameters  used  for 
predictions : 


E  -  35.000  MPa  (Young's  Modulus) 
To  -  1.57 

a  m  3.7  MPa  e’  .  0.055 

S  o 

bx  -  0.35  ba  -  4.0 


o  -  0.25  (Poisson's  ratio) 

Tl  -  0.80 
n  -  0.30 
ba  -  2.0 


Plots  of  mean  pressure  versus  volumetric  strain  In  Figure  5a  show  both 
shear  enhanced  compaction  and  dilatation.  Too  much  dilatation  Is  predicted 
for  the  uniaxial  stress  path  and  the  discrepancy  Is  even  greater  on  the  plot 
of  stress  versus  strain  shown  In  Figure  5b.  This  Is  a  consequence  of  using 
the  measure,  e\  for  all  limit  states.  Predictions  for  the  other  paths  are 

O 

quite  reasonable. 

The  second  example  Involves  foundry  core  which  was  first  proposed  by 
Grady  [22]  as  a  material  for  simulating  cemented  sand  with  material  data 
obtained  by  Waterways  Experimental  Station.  Since  the  density  Is  relatively 
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CL 


Figure 


(a)  Mean  pressure  versus  volumetric  strain. 


.  Experimental  and  theoretical  response  curves  for  a  high-strength 
concrete. 
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Strain,  % 


(b)  Stress  versus  strain. 


Figure  5.  Concluded. 
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low  (1460  kg/m  *) ,  forrestal  proposed  a  modification  to  the  method  of 
sample  preparation  which  raised  the  density  to  1820  kg/m*.  This  density  Is 
similar  to  that  of  dry  porous  tuff  at  the  Sandla  Tonapah  Test  Range  where 
densities  range  from  1620  to  2000  kg/m*.  The  experimental  data  given  here 
were  obtained  by  Terra  Tech  [23]  for  pressures  up  to  400  MPa  on  material 
prepared  at  the  Sandla  National  Laboratory  foundry.  Porosity  was  30  per¬ 
cent.  Material  properties  used  for  predictions  are  given  as  follows: 


E  -  1090  MPa 

yj 

-  0.22 

yo  m  1.06 

\ 

-  0.92 

a  -  1 . 5  MPa 
s 

eo  -  0.55 

n 

»  1.6 

b*  -  0  .IS 

ba  -  0 

b  a 

-  3.5 

The  experimental  results  shown  In  Figure  6  are  composite  curves  from 
two  or  more  tests  with  difference  bars  shown  to  Indicate  the  extent  of  vari¬ 
ation  that  was  observed  among  samples.  Some  of  the  difference  between 
experimental  and  theoretical  results  In  Figure  6a  Is  due  to  the  restriction 
imposed  by  the  functional  form  of  Eq.  (9).  Part  of  the  discrepancies  for 
the  trlaxlal  response  paths  of  Figure  6b  can  be  attributed  to  the  use  of  the 
strain  Invariant  e*  and  the  assumption  that  the  limit  state  is  reached  at 
the  same  value  of  5^  no  matter  what  path  Is  followed.  However  the  results 
are  certainly  adequate  for  engineering  analysis  since  the  variation  In  field 
properties  is  considerably  more  than  that  displayed  by  these  test  samples. 

Overall,  correct  response  features  are  displayed  with  a  rather  elementary 
model  for  an  extremely  compressible  material. 
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5.  CONCLUDING  REMARKS 

This  theory  has  been  developed  under  rather  restrictive  conditions. 

Only  paths  defined  by  monotonlcal ly  Increasing  stress  or  strain  for  positive 
mean  pressure  have  been  considered.  Cyclic  and  rate  effects  have  not  been 
Included  although  It  Is  readily  apparent  that  kinematic  hardening  and  a  rate 
sensitive  flow  surface  can  be  Incorporated  In  the  theory. 

This  Is  one  of  a  large  number  of  Investigations  concerning  the  applica¬ 
tion  of  plasticity  theory  to  frictional  materials.  However,  In  many 
Instances  the  complexity  of  the  models  precludes  their  use  for  analysis  and 
even  for  numerical  computations.  Since  frictional  materials  are  notorious 
for  displaying  experimental  data  with  large  scatter  and  for  being  extremely 
Inhomogeneous  In  the  field.  It  makes  no  sense  from  an  engineering  viewpoint 
to  attempt  to  refine  a  model  Just  to  provide  a  close  match  to  experimental 
data  for  all  stress  paths.  In  the  same  context,  however,  for  general  appli¬ 
cations  It  Is  certainly  necessary  to  have  the  capability  for  predicting  cor¬ 
rect  response  features  which  Include  shear  enhanced  compaction,  dilatation, 
hardening,  and  softening.  The  actual  degree  of  refinement  that  Is  required 
Is  a  matter  of  Judgement  based  on  the  particular  application. 

The  Initial  workers  In  this  field,  as  exemplified  by  the  work  In  Refer¬ 
ences  1  and  2,  were  attacking  essential  problems  with  a  minimum  degree  of 
complexity.  It  Is  In  this  spirit  that  this  model  has  been  developed.  The 
fundamental  concepts  are  not  new.  By  using  stress  Invariants  the  theory  1$ 
completely  general  and  satisfies  the  principle  of  material  frame  Indiffer¬ 
ence.  8y  eliminating  the  second  invariant,  the  theory  Is  simplified  consid¬ 
erably  but  essential  response  characteristics  are  still  predicted.  Funda¬ 
mental  features  such  as  the  limit  state  and  the  Ko  condition  are  displayed 
Inherently  by  the  model.  Compaction  and  dilatation  are  controlled  by  the 
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parameter  y*  In  the  nonassoclated  flow  rule  while  hardening  and  softening 
are  controlled  with  y.  Some  familiarity  with  the  P-L  space  leads  to 
natural  Interpretations  that  are  In  accord  with  experiment. 


To  the  limited  extent  that  the  theory  provides  a  natural  description  of 
physical  phenomena  using  the  clearest  mathematical  structure,  this  paper  Is 
offered  as  an  expression  of  gratitude  to  Professor  E.  F.  Masur,  who  has  pro¬ 
vided  numerous  examples  of  elegant  research  for  others  to  emulate. 
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